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We revisit the use of ideal hydrodynamics to describe bulk- and multistrange-hadron observables 
in nuclear collisions at the Relativistic Heavy Ion Collider. Toward this end we augment the 2+1- 
dimensional code "AZHYDRO" by employing (a) an equation of state based on recent lattice-QCD 
computations matched to a hadron-resonance gas with chemical decoupling at T c h ~160MeV, (b) 
a compact initial density profile, (c) an initial-flow field including azimuthal anisotropies, and (d) a 
sequential kinetic decoupling of bulk (n, K, p) and multistrange (<j>, S, f2) hadrons at T ~ 110 MeV 
and 160 MeV, respectively. We find that this scheme allows for a consistent description of the 
observed chemistry, transverse-momentum spectra and elliptic flow of light and strange hadrons. 
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I. INTRODUCTION 

Experiments at the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC) Q], Q 
suggest that a quark gluon plasma (QGP) is created in 
ultra-relativistic heavy-ion collisions (URHICs) which is 
strongly coupled and behaves like a near-perfect liquid 
[H, H|. In particular, the use of ideal relativistic hy- 
drodynamics enabled a good description of several bulk 
hadron-observables encompassing more than 90% of the 
produced particles [5l-H^. A rapid thermalization of the 
medium, leading to collective phenomena including el- 
liptic flow, could be established and are key to our un- 
derstanding of the macroscopic properties of the fire- 
ball. The success of ideal hydrodynamics and the con- 
clusion that dissipative effects appear to be small has 
more recently led to efforts to quantify these by employ- 
ing second-order viscous hydrodynamics [15l423j . This 
goal clearly requires a good control over any remaining 
uncertainties within the hydrodynamic framework, e.g., 
the equation of state (EoS), initial conditions and imple- 
mentations of the freezeout scenario and/or final-state 
transport. 

Some of these aspects and their interplay, common to 
both ideal and viscous hydrodynamics, are not well un- 
derstood to date. For example, the elliptic flow, v 2 , calcu- 
lated in ideal hydrodynamics seemed to favor an EoS with 
a strong first order phase transition [24[, contradicting 
the finite-T cross-over transition now firmly established 
in lattice quantum chromodynamics (QCD) [2^, [26[. On 
the other hand, the recent progress in solving the so- 
called Hanbury Brown- Twiss (HBT) puzzle required sev- 
eral effects to increase the transverse expansion, includ- 
ing viscosities, initial flow and a hard EoS without phase 
transition [2j . The development of initial flow, prior 
to the thermalization time assumed in hydrod yna mics, 
can be expected on rather general grounds [281 1 2 91 . but 
has only been studied in few works to date [cl l27l [30l — 



|32| |. The initial density profile and initial fluctuations 
are not yet well constrained from first principles, with 
both Glauber- and Color-Glass Condensate (CGC)-based 
approaches currently being investigated [33, HH ■ If dis- 
sipative effects become large, a transition to a transport 
treatment of the bulk is in order, which has been studied 
by coupling hadronic cascades to hydrodynamic evolu- 
tions of the QGP [HI [H, IH, [Hj]. However, it is quite 
possible that the viscosity in the hadronic phase remains 
small for a significant range of temperatures below T c , 
especially if partial chemical equilibrium is implemented. 
The latter becomes problematic in cascade models if the 
inverse of reactions with multi-particle final states need 
to be accounted for, as, e.g., for baryon-antibaryon anni- 
hilation into mesons [36l |37I|. To date, sequential chem- 
ical and thermal freezeouts in URHICs are experimen- 
tally well established, signified by statistical-model fits 
to hadron abundances |3m4(j ] on the one hand (yield- 
ing T c h ps 160 MeV), and empirical blast- wave fits to 
transverse-momentum (pt) spectra of bulk hadrons (it, 
K, p) [H| on the other hand (yielding T fo « 100 MeV). In 
the hadronic EoS figuring into hydrodynamics the num- 
ber conservation of stable hadrons (tv, K, p, p, rj, etc.) 
between T c h and Tf Q can be enforced by introducing per- 
tinent chemical potentials [41] . This has been imple- 
mented into ideal hydrodynamic models 0, 
It was found that bulk-hadron px spectra can be repro- 
duced well at Tf , but the previous agreement with the 
observed elliptic flow, ^(pt), deteriorates [1,|43|, [44|, i.e., 
it is overprcdicted. In viscous hydrodynamics, a system- 
atic investigation of the effects of chemical freezeout on 
bulk observables are still in their beginnings [2l|, HH, H|| . 

The spectra and v-i of multistrange particles have re- 
ceived relatively little attention in hydrodynamic calcu- 
lations thus far. The </>, 5 and ft have no well established 
resonances with bulk hadrons, and elastic i-channel ex- 
change processes are suppressed by the Okubo-Zweig- 
Iizuka (OZI) rule. Therefore, multistrange hadrons are 
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not expected to undergo significant rescattering in the 
hadronic phase and should decouple from the system 
early [13, H3. This is supported by experimental data 
from RHICp], |4iJ . Compared to bulk hadrons mul- 
tistrange particles thus reflect more directly the collec- 
tive dynamics of the partonic stage of the fireball and 
can provide a significant but often neglected constraint 
on hydrodynamic evolution models. For example, in the 
2+1-dimensional hydro-simulations of Ref. [f|, Q~ freeze- 
out has to be carried well into the hadronic phase to be 
compatible with the experimental px spectra. 

In this work we will revisit to what extent ideal hydro- 
dynamics is capable of providing a realistic description 
of the bulk evolution of the medium in Au-Au collisions 
at RHIC. With "realistic" we mean, on the one hand, 
a consistent description of light- and strange-hadron ob- 
servablcs encompassing their abundances, px-spectra and 
V2{pt) at midrapidity for semi-/central collisions. On the 
other hand, we also refer to inputs to, and assumptions 
in, the hydrodynamic treatment which are within the 
uncertainties described above. This includes (i) a state- 
of-the-art equation of state adopted from lattice QCD 
in the QGP phase, matched to a hadron resonance gas 
with hadrochemical freezeout; (ii) a sequential freczeout 
of bulk-hadron chemistry and kinetics, as well as simul- 
taneous kinetic and chemical freezeout of multistrangc 
hadrons (at the universal chemical freezeout); (iii) a com- 
pact initial density profile with relatively large gradients 
and "reasonable" initial radial and elliptic flow. We will 
also elaborate on arguments why the viscosity in the 
hadronic phase of URHICs could be rather small. We 
will, however, neglect the role of initial-state fluctuations, 
which are presumably the main source of higher flow har- 
monics, but do not seem to affect vq.2 much [20| . 

Our computations build on the existing and readily 
available ideal hydrodynamic code package AZHYDRO 
developed by Kolb and Heinz [BJ . One of our goals in per- 
forming a tunc of AZHYDRO is to provide an easily ac- 
cessible yet realistic background medium for quantitative 
applications to electromagnetic and heavy-flavor probes. 
In particular, we do not imply to supersede emerging 
viscous hydrodynamic calculations, which, once fully de- 
veloped, are expected to become the tool of choice. 

Our article is organized as follows. In Sec.|H]we briefly 
recall the basic features of the 2+1 dimensional hydro- 
code AZHYDRO and then focus on its amendments as 
applied in the present work. In Sec. IIII1 we analyze the 
qualitative impact of our amendments on the evolution of 
the bulk medium, i.e., its radial and elliptic flow, and po- 
tential ramifications for hadron spectra, especially freeze- 
out properties. In Sec. lIVl we conduct quantitative fits to 
bulk- and multistrangc-hadron observables, encompass- 
ing multiplicities, pt spectra and elliptic flow. We sum- 
marize and conclude in Sec. |Vl 



II. AZHYDRO AND ITS AMENDMENTS 

The starting point for ideal hydrodynamics (IH) are 
the equations for the conservation of energy and momen- 
tum, 

d„T^ = , (1) 

formulated in terms of the energy-momentum tensor, 
P". As usual, the latter is given by the energy density, 
e, and pressure, p, in the local rest frame assuming ki- 
netic equilibrium, together with a flow field v describing 
the collective motion of the fluid cells. Other conserved 
currents (e.g., for baryon number), can be introduced as 
appropriate, see, e.g., Refs. [flQ for reviews. The system 
has to be closed by specifying an equation of state p(e), 
initialized at a thermalization time To (typically with an 
initial entropy density and flow field) and frozen out in 
the dilute stage (typically at a final energy density us- 
ing the Cooper-Frye prescription [5l[ to convert the fluid 
cells into hadron spectra). In AZHYDRO 0, H| this is 
done in a longitudinally boost-invariant setup leading to 
a 2+1 dimensional evolution. In the following we discuss 
in more detail our modifications to two of the above in- 
gredients, namely the equation of state (Sec. Ill A|) and 
the initial conditions fSec. Ill Bj) . 



A. Equation of State 

The default AZHYDRO code (version v0.2) HQ em- 
ploys an equation of state (labeled "EoS-99" ) consisting 
of an ideal massless quark-gluon gas at high temperatures 
and a hadron resonance gas (HRG) in partial chemical 
equilibrium at low temperatures 1 . The two parts are 
matched via a Maxwell construction at a critical tem- 
perature T c = 165 MeV in a first-order phase transition 
with a mixed phase, see dashed line in the upper panel 
of Fig. [TJ The sound velocity c s vanishes in the mixed 
phase (see the lower panel of Fig.QJ, resulting in a vanish- 
ing acceleration over a relatively long duration. Conse- 
quently, the radial flow of the system largely stalls until 
the end of the mixed phase. In practice, this implies, 
e.g., that multistrangc particles like the Q~ baryon have 
to be decoupled close to the kinetic freezeout of the bulk 
particles, at Tf ~ 100 MeV, to reproduce the pertinent 
experimental spectra at RHIC 

Recent latticc-QCD (1QCD) calculations have now es- 
tablished that, at vanishing baryon chemical potential 
(ps = 0), the transition between the QGP and hadronic 
matter is a smooth crossover [2^, [2(| . The pseudo-critical 
deconfinement transition tem per ature obtained by the 
two leading 1QCD groups (26l. |52| has now converged to 



1 Recently, a patch was released for AZHYDRO by Molnar and 
Huovinen wich also contains a lattice- based EoS 14511. 
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FIG. 1: (Color online) Comparison of pressure (upper panel) 
and speed of sound squared (lower panel) vs. energy density 
as obtained from EoS-99 (dashed lines) and our new equation 
of state (solid lines). 



a common value of around T c ~ 170 MeV [5j|. Some 
differences remain regarding the interaction measure, 
I = e—3p H3, HH • The calculations by the Wuppertal- 
Budapcst (WB) collaboration [55[ reproduce the HRG 
results in chemical equilibrium below T c , while their lat- 
tice action is not ideally suited to the high-T limit where 
their thermodynamic quantities come out slightly below 
the results of the HotQCD collaboration |54[. In our fit 
we decided to employ the WB results in the transition 
region and to smoothly match on to the HotQCD results 
for T > 180 MeV. Once I(T) is specified we calculate the 
EoS, p(e), following the procedure described in Ref. [55| . 
The resulting EoS describes strongly interacting mat- 
ter in thermal and chemical equilibrium. However, in 
URHICs it is well-known that hadron ratios freeze out at 
a temperature of T c h — 160 MeV [38l440| . To account for 
the departure from chemical equilibrium in the hadronic 
phase below this temperature, we follow the approach 
in Ref. [5fjj . by introducing effective chemical potentials 
for hadrons which are stable under strong interactions, 
i.e., pions, kaons, etas, nucleons and antinucleons, in- 



cluding their feeddown contributions (e.g., p — » 2ir with 
Hp = 2/j.tt), usually referred to as "partial chemical equi- 
librium (PCE). The conservation of antibaryon number 
is of particular importance since, in turn, it triggers the 
build-up of pion chemical potentials (56|. We have veri- 
fied that the pertinent chemical-equilibrium EoS is con- 
sistent with the 1QCD results. In the construction of our 
URHIC EoS we then replace the chemical-equilibrium 
part with the PCE part. We will refer to the resulting 
EoS as "latPHG" EoS 2 , which is compared to EoS-99 
in Fig. [TJ As expected, the most notable differences are 
in the transition regime, where the pressure in the lat- 
PHG EoS is enhanced; but, most importantly, the speed 
of sound remains large, c 2 = dp/de ~ 0.15 — 0.20, com- 
pared to zero in EoS-99. As is well known, this will have 
significant ramifications for the radial flow, especially to- 
ward the end of the transition, e ~ 0.5 GeV/fm 3 , and 
in connection with the kinetic freezeout of multistrange 
particles. 



B. Initial Conditions 

Initial conditions currently constitute one of the 
largest uncertainties in hydrodynamic simulations of 
URHICs 0. In default AZHYDRO, a combination of 
wounded-nucleon and binary-collision density is used to 
initialize the entropy density, 

i u\ 4. rn r,r n Bc(x, y;b) n WN (x,y;b) 

s(T ,x,y;b) = const 0.25 n , +0.75 , 

n B c(0, 0;0) n W N(0, 0;0) 

(2) 

at an initial time of t = 0.6 fm/c for each impact pa- 
rameter b. This ansatz leads to good agreement with 
the observed centrality dependence of charged particle 
multiplicities @, Other initial density profiles, e.g., 
inspired by the CGC, have also been used l59j| . 

However, to accelerate the build-up of radial flow, 
which is essential for describing spectra of multistrange 
particles at T c h, a compact initial profile with large ini- 
tial pressure gradients is favored. Such a profile is fur- 
thermore an essential ingredient to a realistic description 
of HBT radii 27J. As a limiting case, we choose the 



entropy-density profile to be solely proportional to the 
binary-collision density 



s(t ,x, y;b) = C{b) n BC (x,y]b) . 



(3) 



While our fits below favor compact profiles they do not 
necessarily dictate collision scaling, albeit it turns out to 
work well. We calculate nee using an optical Glauber 
model 0, H3] • Compact profiles have also been used in 
some other hydrodynamic simulations pll . l43l [60| . Since 



2 Following the systematic analysis of hadron observables in 
Ref. [38| we also introduce a "strangeness suppression" factor 
7s = 7s — 0.85 for each net (anti-) strange quark in a hadron. 
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the particle multiplicity does not scale with iibCi the co ~ 
efficient in Eq. ([3]) needs to become ^-dependent j43| . 

Another uncertainty in the initial conditions concerns 
the possibility of pre-equilibrium flow. In most hydrody- 
namic simulations, the initial transverse collective veloc- 
ity is assumed to be zero. However, it has been argued 
that flow can easily emerge before kinetic equilibrium is 
established [HI, [23] . An initial radial flow field has been 
implemented into a few calculations and shown to im- 
prove the agreement with bulk-hadron data [1, [Hf (e.g., 
by reducing the final V2 in calculations with PCE in the 
hadronic EoS, or by improving on the HBT data). In the 
present work, we go one step further and adopt a non- 
trivial pre-equilibrium flow field including finite clliptic- 
ity. Specifically, we employ the empirical ansatz proposed 
in Ref. [6l|, which was successfully used to fit bulk and 
multistrange observables in a sequential kinetic freezeout 
scenario |62| . The transverse velocity is parameterized in 
terms of the spatial coordinates r and 4>s as 



v(r, 4> s ) = r[a + a 2 cos(2^6)] 



(4) 



where a quantifies the surface radial flow and «2 the 
anisotropy. The azimuthal angle (f>b of the flow vector 
can be tilted away from <fi s through 



R?. 

tan0b = n-j^t&xxps 



(5) 



where k is a tunable parameter dependent on the impact 
parameter b, and i? x = Rq — 6/2, R y = \/ Rq — (b/2) 2 
are the short and long half-axis of the initial ellipse, and 
Ro is the radius of the colliding nuclei (with identical ^4). 
Finally, the normalized radius f in Eq. (j4]) is 



Rl 



(6) 



We keep the initial time ro = 0.6fm/c as in the default 
AZHYDRO. 



III. BULK EVOLUTION 

Before turning to quantitative fits to hadron spectra in 
the following section, we first illustrate the effects of our 
amendments on the bulk-matter evolution of a fireball 
in semicentral Au-Au at RHIC. The bulk evolution in a 
hydrodynamic system can be characterized by the time 
evolution of the average transverse radial flow, (vt), and 
of the anisotropy of the energy- momentum tensor [5| , 



J dxd y (T xx - ryy) 

J dxdy(T™ + Tyy) 



(7) 



In Fig. [3J we compare the longitudinal proper-time evolu- 
tion of both quantities in default AZHYDRO to simula- 
tions with updated EoS and with both updated EoS and 
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FIG. 2: (Color online) (a) Proper-time evolution of the av- 
erage transverse radial velocity, («t), for different hydro sce- 
narios: (i) default AZHYDRO with EoS-99 (green dashed 
curve), (ii) default AZHYDRO with the new EOS (blue dot- 
ted curve), and (iii) the new EoS with initial flow and com- 
pact initial density profile (red solid curve), (b) The same 
comparison for the time evolution of the energy-momentum 
anisotropy, ep. (c) The evolution of the temperature in the 
central cell of the transverse plane. 



modified initial conditions. The parameters for the lat- 
ter correspond to the final data fit discussed in Sec. lIVB] 
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Since the new EoS does not feature a vanishing acceler- 
ation at (or around) T c , the knee in (vt(t)} and the dip 
in £p(r) disappear. However, the acceleration during the 
first 3fm/c is slightly lagging behind in latPHG relative 
to EoS99 due to the non-ideal behavior of the former. 
As a result, the total radial flow is not very different at 
the beginning of the hadronic phase, which also implies 
that the anisotropy keeps increasing significantly in the 
hadronic evolution of both scenarios. These (possibly un- 
desired) properties can be modified by introducing more 
compact initial profiles and initial flow, as shown by the 
pertinent solid lines in Fig. [5J the radial flow increases 
by ca. 50% at T c , while the eccentricity essentially levels 
off thereafter. Both features are crucial in fitting mul- 
tistrange particle spectra at T c , and also improve the 
description of bulk particles at Tf Q (in particular the re- 
duced anisotropy, but the increased final flow also helps 
to describe the pt spectra out to higher px than before). 
The more rapid expansion of the fireball shortens the life- 
time of the fireball by almost 15%, from ~10.6fm/c in 
default AZHYDRO to ~9.0fm/c in the fully amended 
case (both with e do c = 0.1094 GeV/fm 3 ). 

The lower panel in Fig. [2] shows the time evolution 
of the temperature in the central cell of the transverse 
plane. While the initial temperature in the fully amended 
case is higher than in the default AZHYDRO (due to 
more compact initial-entropy profile), it turns out that 
the duration of the hadronic phase in the two cases is 
comparable (note that T c = 170 MeV vs 165 MeV in the 
amended and default AZHYDRO, respectively). 



IV. HADRON OBSERVABLES 

The amended AZHYDRO as described above has been 
utilized to conduct "eye-ball" fits to simultaneously de- 
scribe spectra and i>2 of bulk and multistrangc hadron 
(at Tfo and T c h, respectively) in Au-Au collisions at full 
RHIC energy. We have focused on two centralities, 0- 
5% ("central") and 20-30% ( "semicentral" ) , which, for 
simplicity, we have approximated by fixed impact pa- 
rameters (6=2.3 and 7.38 fm, respectively) corresponding 
to participant numbers, Y par t, calculated in the optical 
Glauber model [38] for the two experimental selections. 
The central initial-entropy density, sg = s(r , 0,0; 6), is 
then adjusted to the experimental hadron multiplicities 
for the two centralities, evaluated at kinetic freezeout in 
the evolution (e fo = 0.1094 GeV/fm 3 , or T fo = 110 MeV) 
incorporating resonance decays (the spectra and feed- 
down from multistrangc hadrons, for which we assume 
early kinetic freezeout, are evaluated at T c h = 160 MeV). 
The two parameters in the initial-flow field basically con- 
trol the radial flow (ao) and, subsequently, the elliptic 
flow (a 2 ). Within our accuracy they turn out centrality- 
independent, provided the parameter for the "tilt" be- 
tween position vector and flow vector is suitably in- 
creased for more peripheral collisions (which is consis- 
tent with intuition, e.g., going to one in the academic 





so (fm 


3 ) To (MeV) 


ao 


012 K 


0-5% 


159.5 


399.9 


0.13 


0.004 1.0 


20-30% 


133.1 


377.9 


0.13 


0.004 4.2 



TABLE I: Initialization parameters for central (0-5%, b = 
2.3 fm) and semicentral (20-30%, b = 7.38 fm) Au-Au 
(y/SNN~200 GeV) collisions; so is the initial-entropy density 
in the center of the transverse plane and To is the perti- 
nent temperature resulting from a collision-density overlap; 
ao, Q!2 and k parameterize the initial anisotropic flow profile, 
cf. Sec.HTBl 



limit of 6=0). The resulting parameter values are sum- 
marized in Tab. HI The kinetic freezeout temperature, 
Tfo = 110 MeV, also turns out to be identical for the 
two centrality classes within our accuracy. 

A few remarks arc in order concerning the applicabil- 
ity of hydrodynamics in the (later stages of the) hadronic 
phase. It is often argued that the latter carries large 
viscosity thus mandating a transport treatment. How- 
ever, there are several arguments suggesting that the 
hadronic evolution in URHICs does not carry large rj/s 
ratios. First, since the QCD transition at /zb=0 is pre- 
sumably close to a second order one, it is likely that 
r]/s possesses a minimum around the pseudocritical tem- 
perature of T c ~ 170 MeV. The question then is how 
fast rj/s increases with decreasing T [63M66| . Dilepton 
measurements (and pertinent calculations) at the SPS 
show that the average p-meson width in the hadronic 
phase is substantially broadened, by more than 200 MeV 
(cf. Ref. 67[ for a recent review). This corresponds to a 
mean-free-path of below 1 fm, and is comparable to the 
thermal kinetic energy of the p, KE ~ 1.5T ~ 225 MeV. 
Even at thermal freezeout the p is still broadened by 
ca. 100 MeV. From another angle, but using similar tech- 
niques (i.e., effective hadronic interactions), charm diffu- 
sion has been evaluated in hadronic matter in Ref. [o8| . 
While the diffusion coefficient (which, in units of the ther- 
mal wavelength, 1/ (2irT), is roughly proportional to T}/ s) 
increases appreciably toward lower temperatures in equi- 
librium matter, the inclusion of effective chemical poten- 
tials only leads to a ^30% increase when going down from 
T=170MeV to 100 MeV. We thus believe that viscosity 
effects in hadronic matter as formed in URHICs may be 
significantly smaller than commonly assumed. 

In the following sections we present the results for 
the observablcs, starting with inclusive yields (hadro- 
chemistry, Sec. I IV A[) and then becoming more differen- 
tial with p t spectra (Sec. IIVBI) and azimuthal dependen- 
cies (Sec- irvCl) . 



A. Particle Yields 

The PCE part of our EoS ensures that the stable- 
hadron numbers are approximately conserved between 
chemical and thermal freezeout. In Tab. [TT] we compare 
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our 7T, K and p multiplicities in central Au-Au colli- 
sions, calculated at T fo , to STAR [H and PHENIX [H] 
data. Since we employ the chemical freezeout temper- 
ature we adjusted the total entropy to obtain the cen- 
tral value of the proton number measured by STAR (in- 
cluding weak feeddown); consequently, the calculated 
7r (corrected for weak feeddown), K and p yields are 
within the experimental errors (residual deviations in- 
side the errors may be caused, e.g., by slight variations 
in the hadronic resonances included in the EoS). Good 
agreement is also found with the kaon measurements of 
PHENIX, while the calculated pion yields are slightly 
above the 1-er upper limit of the PHENIX pions (which 
include weak feeddown). A more significant discrepancy 
arises with the PHENIX anti-/protons (weak feeddown 
corrected). Similar findings are reported in other hydro- 
model fits [ID, [7^ . We note that a slightly higher cen- 
trality selection in the PHENIX data could still be com- 
patible with our calculated kaons while improving the 
agreement with pions and anti-/protons. 



B. Single- Particle Spectra 

Let us first turn to the multistrange hadron spectra (</>, 
S and fl) which we evaluate at the chemical-frcezeout 
temperature, T c h = 160 MeV, where the energy den- 
sity is e c h = 0.372 GeV/fm 3 . As an additional central- 
ity class we consider 20 — 40% central Au-Au collision 
which we approximate with b = 8.04 fm with an initial- 
entropy density in the center of the transverse plane of 
,sq = 130.9 fm -3 to reproduce the pertinent bulk-particle 
yields; all other parameters are the same as for 20 — 30% 
ccntrality (including the strangeness suppression factor 
of 7 S = 0.85 [HI). The comparison of our calculated 
multistrange hadron spectra to STAR measurements are 
shown in Fig.[3j showing good agreement in both central 
and semiccntral collisions. As typical for hydrodynamic 
simulations, the description extends to slightly higher px 
for higher centrality, here up to pt — 4 — 5 GeV for hy- 
pcrons. 

With the same set of parameters, the freezeout of bulk 
particles, (it, K, p) is evaluated at Tf = 110 MeV (at 
an energy density of ef = 0.1094 GeV/fm 3 ). This is 
consistent with values extracted from systematic blast- 
wave fits performed by the experimental collaborations 
to their data [H, [74| • Consequently, our calculated spec- 
tra agree well with their spectra as well, cf. Fig. @] albeit 
with absolute normalization owing to the implementation 
of chemical freezeout. One exception are the PHENIX 
protons (including weak feeddown corrections from hy- 
perons): while the spectral shape is well described, a 
renormalization factor of 0.71 needs to be applied for an 
optimal description of the absolute yields (in accord with 
discussion in Sec. II V Al and Tab. HI)) . 

In Ref. [43[ , which employs an EoS similar to that used 
in AZHYDRO, it was deduced that bulk-particle spectra 
can be fitted better with a chemical freezeout tempera- 



£ 10" 

1 10" 

10" 5 
10"' 




STAR,- ,0-5% " • t 

* STAR,(P_ + +n")/2,0-5% 



,,b=2.3 fm 
H",b=2.3 fm 
fi",b=2.3fm 







2 3 
p T (GeV) 



1 

0.1 

>, 0.01 

-a 

% 1E-3 

Q. 

g 1E-4 
1E-5 



f J 



1E-6 r 



(b) 



■ STAR,<|>,20-30% •* 

STAR,H",20-40% " « 
* STAR,(fi + +fi )/2,20-40% 

<j>,b=7.38 fm I *■ " 

: 3" ,b=8.04 fm 

n",b=8.0.4 fm, . , . , 







1 



2 3 
p T (GeV) 



FIG. 3: (Color online) (a) Our calculated pr-spectra of multi- 
strange particles (<f>, S and fl) in 0-5% central Au-Au collisions 
compared to STAR data [7l|, UM > the thick (thin) lines corre- 
spond to kinetic freezeout at T ch = 160 MeV (T fo = 110 MeV, 
with identical normalization), (b) The same comparison for 
semicentral collisions (without thin lines for late freezeout). 



ture of 150 MeV, ca. 10 MeV lower than the value sug- 
gested by statistical model fits [38M40I ] . Our results indi- 
cate that with a modern lattice EoS and initial flow the 
chemical freezeout temperature in hydrodynamics can be 
made compatible with statistical model fits. 

To illustrate the significance of the early freezeout of 
multistrange hadrons in our calculations, we also plot 
their spectra at Tf Q = 110 MeV in central Au-Au (thin 
lines in the upper panel of Fig. [3]). One sees that the 
additional radial flow developed in the hadronic phase 
hardens the spectra substantially leading to a systematic 
overprcdiction of the data with increasing px- The dis- 
crepancy is largest for the Q ~ and less pronounced for 
the S and 4>. The latter may not be surprising; e.g., 
the S possesses one light valence quark and at least one 
pion-resonance excitation (S(1530)), while the cf> couples 
strongly to both np and KK channels. Thus both <j> and 
5 might develop a reaction rate in hadronic matter which 
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dN/dy 7T+ 7r K + K p p 

STAR 322 ±25 327 ± 25 51.3 ± 6.5 49.5 ± 6.2 34.7 ± 4.4 26.7 ± 3.4 

PHENIX 286.4 ± 24.2 281.8 ± 22.8 48.9 ± 5.2 45.7 ± 5.2 18.4 ± 2.6 13.5 ± 1.8 

Y feeddown included 312.2 314.4 48.2 48.4 34.0 26.2 

Y feeddown subtracted 303.1 303.1 48.1 48.2 25.8 19.9 



TABLE II: Comparison of the calculated bulk-particle yields with STAR and PHENIX measurements at midrapidity in 0-5% 
central Au-Au collisions at ^/snn = 200 GeV. Results are shown with and without feeddown from hyperon decays. The initial 
entropy has been adjusted to reproduce the observed STAR anti-/protons [3^]. The STAR pions are weak feeddown corrected, 
while the protons and anti-protons include hyperon feeddown. PHENIX anti-/protons are weak feeddown corrected [6S|, but 
not the pions. 
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FIG. 4: (Color online) (a) Our calculated pt spectra of bulk 
particles (tt + , K + and p) in 0-5% central Au-Au collisions 
compared to PHENIX [6S| and STAR data (b) The 

same comparison for semicentral collisions. 



FIG. 5: (Color online) (a) Our calculations for the elliptic- 
flow coefficient as a function of pr for <j> and S compared to 
STAR data [fl [7(| in 10-40% central Au-Au collisions; the 
thick (thin) lines correspond to kinetic freezeout at T c h = 
160 MeV (T io = 110 MeV, with identical normalization); the 
calculations employ the same impact parameter as for the 20- 
30% centrality class, (b) The same comparison for tt, K and 
p with STAR data for 20-30% centrality from Ref. [z3]; only 
the results for Tj = 110 MeV are shown. 



allows them to pick up some additional collectivity after 
chemical freezeout, leading to an effective freezeout tem- 
perature slightly below T c h; similar findings are reported 
in Ref. H3 for the 0. 



8 



C. Elliptic Flow 

We finally turn to the elliptic flow, which in our calcu- 
lation is controlled, to a limited extent, by the parame- 
ters ct2 and k characterizing the anisotropy of initial flow, 
and, to a lesser extent, by the interplay with the initial 
radial flow (the initial spatial anisotropy is fixed by our 
assumption of the jiro profile of the entropy density) . As 
is well known 0, @, M [ill, the «2 of massive particles is 
reduced at low pt by a large radial flow. This, in particu- 
lar, also occurs when initial flow fields are introduced Q, 
thus mitigating the problem of a growing pion v-> in the 
hadronic phase when a PCE EoS is employed [8l l43l |44| . 

Our results for the anisotropy coefficient of multi- 
strange and bulk particles in semicentral Au-Au colli- 
sions are compared to STAR data in the upper and lower 
panel of Fig. [5J respectively. Our fit yields fairly good 
agreement for both hadron classes. Clearly, the afore- 
mentioned problem of previous PCE implementations is 
largely resolved due to our more explosive expansion, in 
connection with initial flow fields, which suppress a sig- 
nificant increase of vi in the late stage of the hydrody- 
namic evolution. This was already indicated by the satu- 
ration of the anisotropy of the energy-momentum tensor 
anisotropy in Fig. [2] 

To illustrate the significance of sequential freezeout, we 
display again the v-i of </> and S at Tf Q = 110 MeV by thin 
lines in the upper panel of Fig. [5J Compared to the px 
spectra (the upper panel of Fig. [3]), the current vi data 
are less discriminatory for the freezeout temperature of 
multistrangc hadrons. In accordance with the remarks 
at the end of Sec. IIVB1 the i>2 of S may also favor a 
kinetic- freezeout temperature slightly below T c h. 

V. SUMMARY AND CONCLUSION 

In the present study we have explored the capabil- 
ity of ideal hydrodynamics to simultaneously and quan- 
titatively describe bulk- and multistrangc-hadron spec- 
tra and V2 in Au-Au collisions at RHIC. Specifically, we 
have augmented an existing 2+1D ideal hydro code by 
(i) an equation of state compatible with recent lattice 
QCD data matched to a hadron resonance gas in partial 



chemical equilibrium, (ii) a sequential kinetic freezeout of 
multistrange and bulk particles, and (iii) a compact ini- 
tial density profile with non-zero radial and elliptic flow. 
We deem these amendments "reasonable" in the sense 
of being either suggested by theory (lattice EoS), experi- 
ment (chemical freezeout) and empirical fits (early kinetic 
freezeout of multistrange hadrons), or at least plausible 
(initial flow and compact profiles) . The above items also 
encompass three of the four components identified for 
solving the pion HBT problem. The main practical con- 
sequences of the modifications are a more rapid build-up 
of the radial and elliptic flow, where the latter essentially 
levels off after hadronization. Phenomenologically, the 
underlying parameters (basically the three in the initial- 
flow field parametrization) can be adjusted as to render 
multistrangc hadrons' kinetic freezeout at T c h compati- 
ble with data, and to subsequently avoid an overshoot- 
ing of the bulk-particle t>2 at Tf D . We consider this a 
significant improvement over existing ideal-hydro calcu- 
lations, thereby corroborating the empirical picture of 
early freezeout of multistrangc hadrons. 

We did not include the effects of viscosities nor of 
initial-state fluctuations, which are both clearly needed 
for realistic hydrodynamic simulations of heavy-ion colli- 
sions. However, we believe that our study can still serve 
as a useful baseline for the effects we have included. In 
addition, we anticipate that our amended AZHYDRO 
can be a valuable tool as a realistic fireball background 
for evaluating hcavy-quark(onium) and electromagnetic 
probes. 
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